Optimal control of population transfer in Markovian 
open quantum systems 



Wei Cui*^'^, Zairong Xi*^'*, Yu Pan'^'^ 

°'Key Laboratory of Systems and Control, Academy of Mathematics and Systems Science, 
Chinese Academy of Sciences, Beijing 100190, P. R. China 
^ Craduatc University of Chinese Academy of Sciences, Beijing 100039, P. R. China 



Abstract 

There has long been interest to control the transfer of population between 
specified quantum states. Recent work has optimized the control law for 
closed system population transfer by using a gradient ascent pulse engineer- 
ing algorithm Here, a spin-boson model consisting of two-level atoms 
which interact with the dissipative environment, is investigated. With opti- 
mal control, the quantum system can invert the populations of the quantum 
logic states. The temperature plays an important role in controlling popula- 
tion transfer. At low temperatures the control has active performance, while 
at high temperatures it has less effect. We also analyze the decoherence be- 
havior of open quantum systems with optimal population transfer control, 
and we find that these controls can prolong the coherence time. We hope 
that active optimal control can help quantum solid-state-based engineering. 

Keywords: PACS: 32.80.Qk, 33.20.Bx, 33.80.Be, 42.50.Hz 
1. Introduction 

There exists a widespread belief that future technologies will employ typ- 
ical quantum features of microscopic systems, like entanglement or superpo- 
sitions of quantum states. The most prominent examples are the rapid devel- 
opment of quantum information theory and computation theory , in- 
cluding realizing high-speed quantum computation, and high-security quan- 
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turn communication. Both quantum computation and quantum information 
are the study of information processing tasks that can be accomphshed in 
quantum mechanical systems. In quantum physics, a system composed of 
two or more quantum subsystems can be in a state that has no classical 
counterpart. It is in a superposition of two or more quantum states with 
a precise phase of the whole system. This property is quantum coherence, 
which can show entanglement. Control of quan tum coherence and entangle- 
ment is of great importance 0, 0, 0, SJo, 10, 11, 12, 13|. The control of closed 



quantum systems is well established [16|]. Efforts to extend these studies to 



open quantum systems l3] , where the systems of interest interact with their 
surrounding environments, are now underway. The aim of dynamical con- 
trol in open quantum systems is to suppress effects of the environment in 
order to preserve the quantum properties, including quantum coherence and 
entanglement, etc. The Zeno effect has been used to achieve certain level 
of control on quantum systems, inhibiting the decay of excited states due 
to many repeated measurements. This form of quantum control has been 



studied in several works, including [10|, [18 
Superconducting qubits 



19, 20 



21 



have also been studied as ways to control 
and interact with naturally formed quantum two-level systems in supercon- 
ducting circuits. The two-level systems naturally occurring in Josephson 
junctions constitute a major obstacle for the operation of superconducting 
phase qubits. Since these two-level systems can possess remarkably long de- 
coherence times. References 22|, |23|, |2J] showed that such two-level systems 
can themselves be used as qubits, allowing for a well controlled initialization, 
universal sets of quantum gates, and readout. Thus, a single current-biased 
Josephson junction can be considered as a multi-qubit register. It can be 
coupled to other junctions to allow the application of quantum gates to an 



arbitrary pair of qubits in the system. These results [22|, |23|, |2J] indicate 



an alternative way to control qubits coupled to naturally formed quantum 
two-level systems, for improved superconducting quantum information pro- 
cessing. Indeed, these predictions have been found experimentally in |25| . 
More recently, reference 26| applies quantum control techniques to control 
a large spin chain by only acting on two qubits at one of its ends, thereby 
implementing universal quantum computation by a combination of quantum 
gates on the latter and swap operations across the chain. They 26|] show that 
the control sequences can be computed and implemented efficiently. More- 
over, they discuss the application of these ideas to physical systems such as 
superconducting qubits in which full control of long chains is challenging. 



2 




1. Chirping 




b. Bidirectional transferring 



Figure 1: (Color online)a. Schematic view of inducing population transfer from the ground 
state to the excited state[33|- b. Schematic view of single-qubit logic gate(30|. 



Considerable experimental and theoretical attentions have been paid in 
this field, especially using superconducting quantum circuits. Reference [28 
analyzes the optical selection rules of the microwave-assisted transitions in a 
fiux qubit superconducting quantum circuit. They 28| show that the parities 
of the states relevant to the superconducting phase in the charge qubit can 
be controlled via the external magnetic fiux. For certain values of the fiux, 
the selection rules are the same as the ones for the electric-dipole transitions 
in usual atoms. In other cases, the symmetry of the potential of the artifi- 
cial "atom" is broken, a so-called A-type "cyclic" three-level atom is formed, 
where one- and two-photon processes can coexist. They also study how the 
population of the three states can be selectively transferred by adiabatically 
controlling the electromagnetic field pulses. Different from A-type atoms, 
the adiabatic population transfer in that three-level atom can be controlled 
not only by the amplitudes but also controlled by the phases of the pluses. 
Thus, they achieved a pulse-phase-sensitive adiabatic manipulation of quan- 
tum states in this three-level artificial atom. Usually, only the amplitude was 
considered for adiabatic control. This example of control on quantum cir- 
cuits has been recently studied experimentally [29^. In 30|, a novel approach 
was proposed to coherently transfer populations between selected quantum 
states in one and two qubit systems by using controllable Stark-chirped rapid 
adiabatic passages (SCRAPS) (see Fig.l). These time-insensitive evolution 
transfers, assisted by easily implementable single-qubit phase-shift opera- 
tions, could offer an attractive approach to implement high-fidelity single- 
qubit NOT operations for quantum computing. Specifically, this proposal 
could be conveniently demonstrated by existing Josephson phase qubits. 

However, a quantum system can never be isolated from the surrounding 
environment completely 17|, |3l|. As a result, a randomization of the phase of 
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the quantum system takes place, and the initial quantum state ends up in a 
classical state. Thus, it is an important subject to analyze the quantum de- 
cay induced by the unavoidable interaction with the environment. Population 
transfer control is an important problem with application to various quantum 
systems, and many results are obtained with closed-loop learning control and 
quantum optimal control for large transfer time, respectively. The SCRAP- 



based quantum gates proposed in [3(| evolve in the closed quantum system, 



in which the surrounding environment effect was not considered. Recently, J. 
H. Schonfeldt, J. Twamley and S. Rebic considered the closed system popula- 
tion transfer control in stark- shift- chirped rapid adiabatic-passage (SCRAP) 
technology Their main result is that by the gradient ascent pulse en- 
gineering algorithm the average fidelity of population transfer over a wide 
range of detunings for both the ground to excited state detuning and the 
ground to target state detuning can be improved. In this paper, we make 
use of general optimal control theory to transfer the quantum population for 
quantum computing in Markovian open, dissipative quantum systems. 

The paper is organized as follows. In Sec. II we introduce the con- 
trolled spin-boson model. The Lindblad type master equation for driven 
open quantum systems with the noise and dissipation effects are presented 
in this section. In Sec. Ill, we analyze the optimal control of quantum sys- 
tem dynamics. Both theoretical analysis and numerical demonstration are 
presented in this section. Conclusions and prospective views are given in Sec. 
IV. 

2. Master equation 

A quantum two-level system coupled to an environment is always modeled 
by a spin degree of freedom in a magnetic field coupled linearly to an oscillator 



bath with Hamiltonian 33 



H = Hc + Hb + Hi. (1) 
Here the oscillator bath is described by Hb = Xl/t ^kO^k'^k'-, the spin's observ- 
able az is coupled with the bath "force" operator X = X]fc(5'fc«I + QkCik)- In 
the following, the Planck constant h is assumed to be 1. Let the controlled 
part 

He = -Bz^z + -Bx(^x 

= ^AE (cos r]az + sin r]ax), (2) 
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where the mixing angle t] = arctan(5^/i?2) determines the direction of the 
effective magnetic field in the x — z plane, and the energy splitting between 
the eigenstates is AE = ^ Bl + Bl. B^, B^ are the external controls. For 
the free system Hamiltonian Hq = with cuq = 1 as the norm unit 

throughout this paper, and Hq can be contained in He- Usually, identifying 
the nature of interactions in a quantum system is essential in understanding 
it. Acquiring information on the Hamiltonian can be difficult for many- 
body systems because it generally requires access to all parts of the system. 
Reference [l^l showed that if the coupling topology is known, the Hamiltonian 
identification is indeed possible indirectly even though only a small gateway 
to the system is used. Surprisingly, even a degenerate Hamiltonian can be 
estimated by applying an extra field to the gateway. This information can 
then be used for achieving a better control of the quantum system. 

Under the adiabatic approximation, it is natural to describe the evolution 
of the system in the eigenbasis of He- From Eq.(|2]), the eigenvalues are 
±^AE and the corresponding instantaneous eigenstates are 

|A+(t)) = cos||0) + sin||l), (3) 

|A_(t)) = -sin||0) + cos||l). (4) 

The goal of the theoretical treatment is to analyze decoherence and the 
controlled population transfer, which is given by the elements of the reduced 
density matrix, defined as ps(t) = trB[ptot(t)], where ptot is the total density 
matrix for both the system and the environment, trs the partial trace taken 
over the environment. The effect of the environment on the dynamics of the 
system can be seen as a interplay between the dissipation and fluctuation 
phenomena. 

The simplest quantum system is a two-level system, whose Hilbert space 
is spanned by two states, an excited state |e) and a ground state \g). The 
Hilbert space of such a system is equivalent to that of a spin-| system. 
The corresponding Pauli operators are cxi = \e){g\ + \g){e\,a2 = —i\e){g\ + 
i\g){e\,a3 = |e)(e| — \g){g\, satisfying the commutation relations [ui, cjj] = 
2ieijkO'k, and the anticommutation relations {ai,aj} = 26ij. Then the con- 
trolled master equation of two-level system has the following Lindblad form 
[l7t (p instead of ps for simplicity), 

dp{t) i^7o(iV + l)r + + 
= -t[Hc,p] + <^ 2ap(7+ - 
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1 loN f 1 

—pa^a >H —< 2a^pa — aa^p — paa^ > , (5) 

where N = l/(e^sT _ i) denotes the mean number of quanta in a mode with 
frequency uq of the thermal reservoir, 70 is the spontaneous emission rate. 
The matrix elements poo = Pg(t) and pu = Pe(t) are the population of the 
ground and excited state levels, respectively. The off-diagonals poi(^) = Pio(^) 
are the coherences. Let x{t) = Tr[(Tp{t)], which is called as Bloch vector and 
defined by 

Xlit) = poi{t) + pio{t), 

X2{t) = i{pio{t) - polit)), (6) 
Mt) = Pooit) - pn{t). 

Then the Bloch vector form of (I5l) is 

Mt) = -fi.Xi(t)-^70X2(t) + S,X3(t), (7) 

X3(t) = -5,.X2(t)-(2Ar + l)7oX3(t)-7o, 



which can be written compactly as 



where 



and 



x{t) = A{t)x{t) + B{t) (8) 

-^70 5. 
A{t) = I -B, -^70 

-B^ _(2iV + l)7o 

B{t) = 




3. Optimal control of population transfer 

3.1. Optimal control formalism 

Recent work [Ij] has optimized the control law for closed system popula- 
tion transfer by using a gradient ascent pulse engineering algorithm. They 
shows that the optimized pulses perform at a higher fidelity than the stan- 
dard Gaussian pulses for a wide range of detunings (i.e., large inhomogeneous 
broadening). The theory of optimal control was introduced in the theory of 
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34 
35 



automatic control in the 1960s for electrical engineering applications. Bell- 
man and Pontryagin, the two famous scientists, paved the way in this field, 
resp ectively. Prof. Belavkin introduced it to quantum mechanical in 1983 
In quantum chemical, quantum optimal control made great success 
|38| . In the optimal control framework, one starts by defining a 
cost functional which has the function of optimality criteria. Moreover, the 
cost functional will vary from one experimental trial to another, and must 
be thought of as a random variable depending on the measurement output. 
The strategy is then to minimize this cost functional while satisfying the 
constraints of the underlying dynamic equations governing the evolution of 
quantum states, e.g., the master equation and the initial state condition. 
The calculation of the necessary optimality conditions for this optimization 
problem results in a system of coupled equations to be solved. For detail one 
can refer to the papers j^, 0, [l3, 15|. Obviously, the evolution of the state 
variable x{t) governed by the master equation ([7]) depends not only on the 
initial state Xq but also on the choice of the time-dependent control variable 
u{t) = (B^it), Bzit))^ . In this section, we are going to control population 
transfer and suppress the unexpected effect of decoherence by optimal con- 
trol technique that wants to force the system evolving along some prescribed 
cohering trajectories. For some target state a;°(t), let the cost functional as 

J[u{t)] = ^[x{t^),x'^{ts)] + r L{x{t),x''{tf),u{t))dt, (9) 



where the functional represents some distance between the 

system and objects at final time and the functional J^^-^ L{x{t)^x^{tf),u{t)) ac- 
counts for the transient response with L{x(t) , x^ (t j) , u(t)) > 0. The optimal 
control problem considered in this paper is to minimize the cost functional 
J[u{t)] with dynamical constraints ([7]) and initial state constraint a;(0). Using 
the Pontryagin's maximum principle, the optimal solution to this problem 
is characterized by the so-called Hamilton- Jacobi-Bellman(HJB) equation. 
We use this method to solve the population transfer and decoherence in the 
following. 
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Figure 2: (Color online)Tinie evolution of occupation probabilities for N e [0, 10] without 
control. 



3.2. Evolution without control 

At first we consider the free evolution of tlie system ([5]), where the external 
control field = 0, = ooq- Then the free evolution is 

Xi(t) = -^-foXi{t)+UJoX2{t), 

X2{t) = -UoXi{t)-^JoX2{t), (10) 



xs{t) = -(2Ar+l)7oa;3(t)-7o, 



whose solution is 



Xi(t) = e 


- 2^ ^°*[X2(0) sin(wo^) + xi(0) cos(wot)], 


(11) 


X2{t) = e~ 


' 2"" ^''*[X2(0) cos(u;ot) - a;i(0) sin(u;ot)], 


(12) 


xs{t) = e 




(13) 



We observe that the population-component x^ of the Bloch vector decays 
exponentially with rate —{2N + 1)70, while the coherence Xi_2 decay with 
rate —{2N + 1)70/2. The stationary solution is 

xl =x'2 = 0, x' = , (14) 

12 5 3 2A^+1 

and the populations of the lower and upper level are found to be 

p,{t) = pooit) = ^(1 + X3(t)) (15) 



ig-(2Af+l)7ot 



^ 2poo(0) - 1 



2A^ + 1 
8 



1^ 1 



2 V 2A^ + 1 



N=0.01 




5 10 15 20 25 



Figure 3: (Color online)Tinie evolution of occupation probabilities poo(red) and pii(blue), 
with system evolution without control (the dashed line) and control target trajectory (the 
dotted line). 



p,(t)=pn(t) = ^(1-X3(t)) (16) 



_ig-(2Af+l)70i 



2N +1 



+ 2pii(0)-l 



1 



1 



2 V 2N 



with the stationary populations |(1 ± 2aF+t)' respectively. As fc^T — )■ the 

mean number of quanta = l/(e^s^ — 1) — )■ 0, which means that poo(^ ~^ 
oo) — )■ and puit — )• oo) — )■ 1. This is a population transfer process. The 
elementary logic gates in quantum computing networks can be implemented 
by this transfer. However, when > poo(^ oo) = |(1 — 2n+i ^ ^ 
Pii{t — ^ oo) = |(1 + 2N+i ) < Furthermore, as ksT — )• oo the stationary 
populations are Pqq{N oo) = pfi{N oo) = |, which fails to transfer the 
populations. The evolution of populations poo = |(l + ^3(^)) (left) and pn = 
i(l — X3(t)) (right) is plotted for A^ G [0, 10] in Fig. 2, where the dissipation 

constant 70 is set as 0.1, and the the initial state x(0) = (^, 1). In the 
following, choosing the trajectory with A^ = as the target we will use the 
optimal control technology to make the single qubit rotation that completely 
inverts the populations of the quantum states. 
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N=0.2 




Figure 4: (Color online)Tinie evolution of occupation probabilities poo(red) and pii(blue), 
with system evolution without control (the dashed line) and control target trajectory (the 
dotted line). 



N=o.5 




Figure 5: (Color online)Tinie evolution of occupation probabilities poo(red) and pii(blue), 
with system evolution without control (the dashed line) and control target trajectory (the 
dotted line). 
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3. 3. Controlled population transfer 

Precise design of switching on/off logic gates of quantum systems is a 
main problem in quantum computer research. In the following we consider 
the optimal control of population transfer in a dissipative open quantum 
system. For simplicity let the cost functional be: 

J[u{t)] = r \{x{t) - x%t)f + eu^(t)u{t)]dt (17) 

J to 

where 6' > is a weighting factor used to achieve a balance between the 
tracking precision and the control constraints. u{t) = {Bxit), Bzit))'^ is the 
external control. x^{t) is the target trajectory 

/ e-^*(x2(0)sin(5,t) + xi(0)cos(5,t)) \ 
x\t)=\ e-^*(x2(0)cos(S,t) -a;i(0)sin(5,t)) , 
V e-T»*(x3(0) + 1) - 1 / 

which is the solution of the master equation ([5]) in the case of iV = without 
control, i.e. ^ = —il^aztp] + '^{2apa~^ — a^ap — pa~^a}. The dissipator 
of the equation describes spontaneous emission process (rate 70) without 
thermally induced emission and absorption process (A^ = 0, or T = 0). That 
was the reason we chose it as the target trajectory. We will choose it as the 
tracking object in both optimal population control in the following. 
The corresponding Hamiltonian function is 

'Hix{t),u{t),^{t),t) 
= [(x(t) - x\t)f + eu^it)u{t)] + A(t)^[A(t)x(t) + B{t)] 
= [(xi(t) - x%t)y + {x,{t) - xl{t)f + {x,{t) - xl{t)f + e{Bl{t) + Blm 
2N + 1 2N + 1 

+ ^lif)[ ^ 70X1 (t) + B,X2{t)] + \2{t)[-B,Xi{t) ^ 7oX2(t) 

+B,x^{t)\ + A3 (t) [-5,0:2 (t) - (2iV + l)^ox^{t) - 70], 

where \{t) = {Xi{t) , X2{t) , X^it))-'^ is the so-called Lagrange multiplier. The 
optimal solution can be solved by the following differential equation with 
two-sided boundary values, 

= ^ = MtMt) + B{t), 

m = -f = -2[x(t)-x°(t)]-A(t)^A(t), . . 

x(0) = xo, ^ ^ 
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Figure 6: (Color online)Tinie evolution of occupation probabilities poo(red) and pii(blue), 
with system evolution without control (the dashed line) and control target trajectory (the 
dotted line). 



Together with 



we have 



Ai 

A2 

As 



-^702^1 + B,X2, 

-B^X2 - (2A^ + 1)702:3 - 7o, 
-2(xi - x? 

— 2(^X2 — ^2 

-2(x3 - xl 

x(0) = Xq, 



+ ^7oAi + 5.A2, 
- 5,Ai + ^7oA2 + SxA 
-5,A2 + (2iV+l)7oA3, 
\{tf) = 0. 



3; 



m 

du 



dn{x*{t),u*{t),i){t),t) 

du 



(19) 



(20) 



(21) 



Bx{t) = jg{hX2 - >^2X3}, 

= ^{A22;i - Aia;2}- 

In general, no analytic solution of the complicated nonlinear equations exists. 
The numerical demonstration to this problem is considered in the following. 



3.4- Numerical demonstration and discussions 

As example we show in Fig. (3-8) the populations as a function of time 
for the open quantum system with six different values of iV(temperature). In 
our simulations, the system parameters are chosen as: a;(0) = {\/2, y/2,1), 
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Figure 7: (Color online)Tinie evolution of occupation probabilities /3oo(red) and pii(blue), 
with system evolution without control (the dashed line) and control target trajectory (the 
dotted line). 



N=10 




5 10 15 20 25 

Figure 8: (Color online)Tinie evolution of occupation probabilities poo(red) and pii(blue), 
with system evolution without control (the dashed line) and control target trajectory (the 
dotted line). 
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dissipation constant 70 = 0.1, and system frequency cuq = 1 as the norm unit. 
We consider the two-mode control, physically, -B^(t) time-dependent Stark 
shift and Bx(t) external resonant control field acting on the open system. 
In our simulations, the time evolutions of the lower level population Pg{t) = 
Poo(^) = |(1 + 2^3(^)) colored red and the time evolutions of the upper 
level population pg{t) = pn{t) = |(1 — X3{t)) are colored blue, and the free 
evolutions are plotted with dashed line, the control target trajectory with 
dotted line, optimal populations control with solid line. 

More precisely, we study the temperature as a key factor in popula- 

tion transfer control. In fact, N = l/(e'=B^ — 1), which means that T — 
hiuo/{kBln{l + j^)), where ks = 1.380662 x lO-^V • K'^ is the Boltz- 
mann constant and h — 1.0545887 x 10~^^J • s is the reduced Planck con- 
stant. When is from to 10, the temperature changes significantly, 
(0 ~ 8.0182) X lO^^^Wo (K). The smaller the N the lower the tempera- 
ture, the larger the N the higher the temperature. Looking at Figs. (3- 
8)(where we set N by 0.01, 0.2, 0.5, 1, 2, 10, respectively), we can approx- 
imately say that when the temperature is low, for example, N < 0.5 (T = 
0.6956 X 10~^^ujq (K)), we can achieve very good population transfer fidelity. 
And the lower the better. Increasing A^, smaller population transfer were 
observed in Fig. (6) (A^ = 1) and Fig. (7) (A^ = 2). This is because in higher 
temperature regime the energy level spacing is decreased, which results in 
low transfer rate. When A^ = 10, the population is wild oscillations making 
the optimal control failure. It means that the environment induced fluctua- 
tions is large enough, so that the control field is negligible comparing with 
the high-frequency harmonic oscillators of the reservoir. From Figs. (3-8) we 
can sum up three rules, (i) The external control field played an important 
role in population transfer for quantum computing, (ii) The smaller the N 
(temperature) the better the control performance, especially when A^ < 0.2 
the completely population transfer can be achieved, (iii) When temperature 
is high enough the transfer is uncontrollable. 

3.5. Decoherence behavior 

In the following we make efforts to study the decoherence behavior of 
the open quantum systems with optimal population transfer control, briefly. 
The persistence of quantum coherence is relied on in quantum computer, 
quantum cryptography, quantum teleportation, and it is also fundamental in 
understanding the quantum world for the interpretation that the emergence 
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Decoherence without control 




20 



Figure 9: (Color online)Time evolution of decoherence for N e [0, 10] without control. 



of the classical world from the quantum world can be seen as a decoherence 
process due to the interaction between system and environment. From (fTOjl . 
one can easily conclude the decoherence factor 



1 

? _ 2N+1 
2^ ' 



xl{t)+xl{t) 

^^otA(O), 



(22) 



decaying exponentially to with rate — 2^7o- In Figs. (9, 10) we plot the time 
evolution of decoherence for continuous variation of (temperature). Fig. 
(10) shows that at low temperature the decoherence can be well controlled, 
i.e., the decoherence time can be delayed and its amplitude can be amplified. 
With increasing (the temperature), the environment induced fluctuations 
will be large enough to neglecte the control field. Comparing with Fig. (9) 
we find that these optimal controls can also prolong the decoherence time 
and satisfy the controlled quantum logic gate operations, especially in low 
temperature regime. 



4. Conclusions 

In the present work, we have studied the controlled population transfer 
between selected quantum states in the Markovian open quantum system, 
which can generate universal logic gates for quantum computing. We use the 
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Optimal decoherence control 




Figure 10: (Color online)Tinie evolution of decoherence for N £ [0, 10] with optimal 
control. 

optimal control method to control the population transfer. Our numerical 
results indicate that the occupation dynamics behaves differently for the dif- 
ferent environmental condition. The result can be summed up in three rules: 

(i) The external control field plays an important role in population transfer. 

(ii) The smaller the N (temperature) the better the control performance, 
especially when N < 0.2 the completely population transfer can be achieved. 

(iii) When temperature is high enough the transfer is uncontrollable. Solid 
state qubits offer remarkable advantages due to their scalability and control- 
lability. Therefore, optimal control of population transfer perhaps provides 
an attractive approach to generate universal logic gates. We hope such op- 
timal control techniques to experimentally implement quantum elementary 
logic gates in quantum computing networks in the near future. 

In this paper, we restrict the discussions to Markovian open quantum 
systems and show the validity of our optimal control strategy to the popu- 
lation transfer, and we analyze corresponding decoherence behavior and we 
find that these controls can prolong the decoherence time and satisfy the con- 
trolled quantum logic gate operation in low temperature regime. To make 
the model more realistic in the experiments and reduce the low temperature 
constraints it is worth extending the scope to non-Markovian open quantum 
systems where the seemingly lost information can return to the system at a 
later time, which is our further work. 
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